Genuinely Multipartite Concurrence of iV-qubit X-matrices 
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We find an algebraic formula for the iV-partite concurrence of N qubits in an X-matrix. X- 
matricies are density matrices whose only non-zero elements are diagonal or anti-diagonal when 
written in an orthonormal basis. We use our formula to study the dynamics of the iV-partite 
entanglement of N remote qubits in generalized A-party Greenberger-Horne-Zeilinger (GHZ) states. 
We study the case when each qubit interacts with a partner harmonic oscillator. It is shown that only 
one type of GHZ state is prone to entanglement sudden death; for the rest, A-partite entanglement 
dies out momentarily. Algebraic formulas for the entanglement dynamics are given in both cases. 



I. INTRODUCTION 

Even though entanglement was promoted by 
Schrodinger [l| as a fundamental aspect of quan- 
tum theory, and in mathematics it predates quantum 
mechanics by decades 0, its value as a resource for 
a wide range of potential applications was not appre- 
ciated until recently [H, Furthermore, although its 
importance is now largely recognized, witnessing and 
quantifying the entanglement of arbitrary mixed states 
are still open questions. In only few-party cases do 
prescriptions exist for determining the entanglement of 
a mixed state 0-0] • 

The problem becomes much more difficult for multi- 
partite entanglement, entanglement shared between more 
than two parties. Multipartite entanglement is thought 
to play an essential role in many phenomena including 
quantum metrology Q and quantum phase transitions 
[9( . Previous studies of the dynamics of multipartite en- 
tanglement have utilized measures that fail to capture ex- 
actly when multipartite entanglement disappears fl0l - [l3| . 
Additionally there have been many attempts to solve the 
problem of determining the multipartite entanglement of 
a given state fl4Hl6| (see also references in the paper 
by Ma et al. [17]). but an algebraic and/or numerically 
efficient prescription has not yet emerged. An algebraic 
prescription would be especially desirable since it can po- 
tentially open the door for a wide range of analytical in- 
vestigations of entanglement dynamics. 

Recently, based on previous works by Pope and Mil- 
burn [l8j and Love et al. fl9j , a new measure of multi- 
partite entanglement, called genuinely multipartite (GM) 
concurrence, has been proposed, and it has been shown 
that GM concurrence is an entanglement monotone fl7"j . 
This measure reduces to the I-concurrence [2(| for a bi- 
partite system, and to Wootters's original concurrence 
[5| for two qubits. Additionally, an algebraic formula for 
a lower bound of the GM concurrence has been found 
by Ma et al. [l?}- The lower bound, when calculated for 



a two-qubit X-matrix, matches the value of Wootters's 
concurrence. 

The X-matrix of Yu and Eberly [2l[ is a density ma- 
trix of N qubits, written in an orthonormal product ba- 
sis, whose non-zero elements are either diagonal or anti- 
diagonal. The concurrence of a two-qubit X-matrix takes 
a very simple form pH ] and that is why these two-qubit 
states have been extensively used in studying the dynam- 
ics of entanglement between two qubits in many scenarios 
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Restricted forms of X-matrices of more than two qubits 
have been used in some recent studies of the dynam- 
ics of multipartite entanglement. The entanglement mea- 
sures utilized in these studies yield qualitative informa- 
tion about the multipartite entanglement (2f| Hfil • Only 
for Greenberger-Horne-Zeilinger (GHZ) states that un- 
dergo pure dephasing, has there been a successful at- 
tempt that uses a geometric measure to give the exact 
dynamics of A^-partite entanglement [27j ■ 

In view of the fact that the GM concurrence lower 
bound, derived in Ref. [l7[, matches the exact value of 
concurrence for a two-qubit X-matrix, one might wonder 
if the lower bound might also be exact for more than a 
two-qubit X-matrix. In this paper we will prove that this 
conjecture is correct. The lower bound provided by Ma 
et al. [l?} is realized by X-matrices. We thus present an 
algebraic formula for the GM concurrence of an A^-qubit 
X-matrix. This is our principal result and it enables ana- 
lytic formulation of dynamics of A^-partitc entanglement 
in different scenarios. 

We will use our result to resolve the dynamics of 
A^-partite entanglement of Af-qubit GHZ states, which 
are superpositions of two bit-flipped basis states, when 
each qubit (two-level atom) interacts locally with an 
additional system. We assume the additional systems 
are single-mode cavities, initially in the vacuum state 
[iHHU. We will identify two kinds of GHZ states, based 
on whether or not they experience entanglement sudden 
death (ESD) [Hj]. We will show that most of the GHZ 
states do not experience ESD and only the GHZ state 
for which, initially, the two-level atoms are all excited or 
all in their ground states experiences ESD. In each case, 
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we will present algebraic formulas giving the evolution of 
the GM concurrence. 



II. GENUINELY MULTIPARTITE 
CONCURRENCE 

In this section we give a brief introduction to TV-partite 
entanglement and GM concurrence. 

TV-partite entanglement is defined by its opposite, 
biseparability. A pure TV-partite system \ip) £ Hi 65 
H2 ■ ■ ■ <S> T~Ln is biscparablc if there is a bipartition of the 
TV parties Hi ® H2 • • • ® Un = Ha ®Hb, where Ha = 



H n ®H 3 



Hi 



)Hi 



H jN for 



which \ip) = \iPa) <& \iPb), and |^A> G and \tp B ) G H B 
[HUHl- In other words, a pure state is biseparable if it has 
at least one pure marginal (reduced density matrix). An 
TV-partite state that cannot be written as an ensemble of 
biseparable states is an TV-partite entangled state. 

Before introducing GM concurrence, let us introduce 
the set of all bipartitions of TV parties. Each biparti- 
tion is a division of the set {1,2,..., TV} into two non- 
overlapping and non-empty subsets. We call the set of all 
such bipartitions J = {Ji, J2, . . . J2«-i_i}. For example 
for TV = 3 



Ji = {l|2,3}, J 2 
J = {Ji, J2, J%}- 



{2|1,3}, J 3 = {3|2,1}, 



(1) 



Now if we have a pure state \ip), to each of the ele- 
ments of J we can associate two reduced density matri- 
ces A(\tp), Jj), and B(\ip), Jj) by tracing out either of the 
subsystems associated with that bipartition. The bisep- 
arability of a pure state can be determined by whether 
for any of the elements of J, A(\i/}),Jj) and B(\ip),Jj) 
are pure. If so, \ip) is biscparablc. Thus the purity of 
the j th bipartition, to be denoted by Uj(\ijj)), is a key 
parameter. 

For a pure state genuinely multipartite concurrence 
(GM concurrence) is then defined [17| as below. 



C GM (\i')) ■■ = min^l-n.dV)) 

Clearly, Cgm(IV')) > and it is equal to zero if and only 
if \ip) is a biseparable state. For a bipartite system this 
definition reduces to the I-concurrence ■ 

To determine whether a mixed state p is biseparable 
or not one has to determine whether p can be written 
as a convex sum of pure biseparable states. Thus one 
has to check all the ways p can be written as a convex 
sum of pure states (all pure state decompositions). Let 
us distinguish different pure state decompositions of p by 
assigning a superscript, a, to each decomposition 



P 



(2) 



To determine whether a particular pure state decompo- 
sition is a sum of biseparable states or not, we can cal- 
culate the average pure state GM concurrence for that 



particular a 

c a {p) = Y,PiC GM m)) 



= £ Pf I min v^/l-II^k/?)) 



(3) 



Now we are ready to extend the definition of GM con- 
currence to all mixed states: 



Cgm(p) = min C a {p). 



(4) 



If Cgm(p) = this means that there is an a for which 
C a {p) = 0. Then p can be written as a sum of pure 
biseparable states, so p is biseparable. If Cgm(p) > 0, 
there is no a for which the )'s are all biseparable and 
thus p is an TV-partite entangled state. The GM concur- 
rence of TV parties, as defined in Ref. [17j , is a monotone 
of genuinely multipartite entanglement; it distinguishes 
between biseparable and TV-partite entangled states, is 
convex, invariant under local unitary transformations, 
and non-increasing under local operations and classical 
communications (LOCC) [l7| . 

In [ijj, Ma et al. presented a lower bound for the GM 
concurrence. Now, we examine this lower bound for an 
X-matrix of TV qubits. If the orthonormal basis for the 
X-matrix is {|0, 0, . . . , 0), |0, 0, . . . , 1), . . . , |1, 1, . . . , 1)}, 
then one can always write an X-matrix in the form given 
below 



X 



a 2 



-2 



- // 



V 



(5) 



bj 



where n = 2 1 , and we require \zi\ < \/aibi and X)i( a i + 
bi) = 1 to ensure that X is positive and normalized. 
One can see that the matrix looks like an X, explaining 
why density matrices in this class are called X-matrices. 
Without loss of generality we can assume that \fa\b\ > 
■Jaibi. The lower bound of GM concurrence, derived in 
[l7| , for an X-matrix is given by 



C GM > 2max{0,| 



0,1,.. 



(6) 



where Wi = Y^j^i V a j bj ■ In the next section wc will 
show that this lower bound is exact for all X-matrices. 



III. CONCURRENCE OF TV-QUBIT X-STATES 

In this section wc find an algebraic formula for the GM 
concurrence of an N-qubit X-matrix. 
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Before going forward let us introduce some notation. 
Since we are working with qubits, we can represent each 
vector (ket) of the above basis as a number from 
to 2 N — 1 written in the binary basis. For example, 
[0,0, - - - ,0) = |0>, |0, ■ -- ,0,1) = |1>, |0, -•• ,1,0) = |2), 
and so on |1, • • • , 1, 1) = \2 N — 1). We also denote the 
bit-flipped states in the same way, |i) = \2 N — i — 1), 
e.g., |0) = \2 N — 1). In places where we need to label the 
individual qubits, we will do so by using a subscript on 
the bits, 

|0i,0 2 ,0 3 ,l 5 ) = |0>i®|0>a®|0> 3 ®|l>6. (7) 



Since we have assumed that VOi&i > \J a^, it is easy 
to show that \zi \ — Wi < for i > 1, so that Eq. [5] reduces 
to 

C GM (X)>2max{0,\z 1 \-w 1 } (8) 

We will show that this bound is actually an equality. 
First let us prove a lemma that we will utilize in our 
proof. 

Lemma 1. The GM concurrence of an X-matrix for 
which a±bi > aibi, and aj = bj = for all j ^ {i, 1}, 
is 

C GM (X li )=2ma,x{0,\z 1 \-\^ib~} (9) 

Proof. We already know that this quantity is a lower 
bound of GM concurrence. Thus we only need to show 
that it is also an upper bound. We will do this by map- 
ping X\i to a two-qubit density matrix, R, and then show 
that Cgm{Xh) is bounded from above by Wootters's 
concurrence of R, where C(R) = 2max{0, \z\ \ — y/aibi}. 

We perform the mapping by focusing on a specific bi- 
partition of the qubits. The four non-zero diagonal ele- 
ments of Xu are {ai, Oj, &i, bi}, corresponding to projec- 
tors {|0)(0|, |0)(0|, and \1~T)(1~T\} respec- 
tively. Those qubits that contribute 1 to the ket \i — 1) we 
designate as party F. The rest of the qubits we denote 
as party G. For example, with 7 qubits, which we denote 
as {lii 92, 93) 94) 95) 96) 97), where i = 6, the basis states 
are 

|0> =|0 1 ,0 2 ,0 3 ,0 4 ,0 5 ,0 6 ,0 7 ), 
|5) =|0i,0 2 ,0 3 ,0 4 ,l 5 ,0 6 ,l 7 ), 

|127) =|ll) I2,l3,l4,l5, 16, 17>, 

|122) =|li,l2,l3,l4,0 5 ,l 6 ,0 7 ). (10) 

Then party F is given by qubits (q±, g 2 , 93, 94, 96), an d the 
remaining two qubits, (95,97), make party G. Under this 
bipartition, we can write Xu using the following basis 
states, 

I If) = |0i,02,0 3 ,04,0 6 ), I Tf) = |li, 1 2 , 1 3 , 1 4 , 1 6 ), 

Ug) = {|0 5 ,0 7 ), |tG> = |l6,l7>, 



X u = ai| |_f|g)(|fIg I +h | ti^tcXtFtc I 

+ Oi\ lFt G )(lFt G \+bi\ ti4 G )(ti4G I 

+ zi\ IfIgXTfTg I +z*\ TfTcXIfIg I 

+ Zi\ IfIgWfIg I +<l TfIgXIfTg I, (11) 

We see that if we restrict attention to the subspace de- 
fined by the non-zero elements of Xu we can map Xu 
to a two qubit density matrix, R, which, in the basis 

{| IfIg), I 4-fTg) 

,| Tf4-g),| ti^tc)}, reads 

/ ffli Z\ \ 

Xu^R=\ % I ■ (12) 



Now that we have a two-qubit density matrix, we can 
take advantage of Wootters's concurrence. Note that 
from each pure-state decomposition (PSD) of R one can 
make a PSD of Xu by mapping the basis states of the 
two-qubit system back to the multi-qubit basis states. 
We pick the PSD whose average concurrence is the min- 
imum amongst all possible PSD's of R. Thus 

£ = X>hWil, c(h) = X>c(|V«». (13) 

i i 

In Eq. [T31 C(R) is Wootters's concurrence, which, by 
definition, is equal to the minimum average concurrence 
over all possible PSD's of R. As mentioned before, each 
pure state can be mapped back to an N-qubit state 
(\t(>i) -> 1*4), producing a PSD for X u , 

X M = 5>i|*i)<*i|. (14) 

i 

Since GM concurrence is convex by definition, we have 
C G m{Xu)<Y.P^ C gm{\^ 1 ))- (15) 

i 

For a pure state the GM concurrence is defined by 

C GM (\Vi)) - minV^l- 11,(1*,)), (16) 

where the minimum is taken over all bipartitions, J, of 
the iV qubits. Therefore, the GM concurrence must be 
bounded by any specific bipartition, including the bipar- 
tition of the TV qubits to party F and party G. 

c GM {\*i)) < v^/i-n F | G (|*i)). (17) 

Using the same mapping as done for Eq. [T2] it is easy 
to show that C(\ijji)) is equal to the right hand side of 
Eq. [17l Therefore we conclude that 

C G m(Xu) < X>Cgm(|*<» < Y,PiC{\A)) 

i i 

= C(R) = 2 max{0, \z x \ - y/a*bi}, (18) 



4 



where the right most equality is found by evaluat- 
ing Wootters's concurrence for R under the assumption 
a±bi > ciibi. This upper bound matches the lower bound 
and therefore it is the exact value of Cgm(Xu). □ 

Next, we generalize this result to all the X-matrices. 
We do so by decomposing the X-matrix into a convex 
sum of Xu matrices. Let us first look at the case for 
which \zi\ — Wi > 0. 

(a) - wi > 0. Note that \z,\ - Wi < for i > 2 
since \faib~i > \/aibi > \zi\. First by a change of phase 
of the basis, which is a local unitary transformation, we 
change z\ to \z\\. This only changes the phase of the 
other off-diagonal elements. Then we decompose X in 
the following form. 



X = A 



(19) 



i>l 



where Si is an Xu matrix whose two-qubit counterpart 
reads 



Ri 




K bi 




(20) 



where 



An =ai(l- 



\/aib- 



:)• 



A 2 



Vaih 

M,2n = A 2n ,l = \zi\ - Wi, 

Aij = i^{l,2n},or j + {l,2n}, 
ai y 1 (Libi 




(21) 

(22) 

(23) 
(24) 

(25) 
(26) 



It can be shown that S^s are all proportional to valid 
density matrices, since they are non-negative hermitian 
matrices. The proportionality constant is between zero 
and one and can be interpreted as probability. Us- 
ing Lemma [T] one can show that all Si's are bisepara- 
ble matrices (though not normalized). Regarding the 
first matrix in the decomposition, the proportionality 
constant is An + A2 n ,2n 1 and its GM concurrence is 
2(|zi| — wi)/(An + A 2 n : 2n)- Due to the convexity of 
GM concurrence we conclude that 



A 2 



2(|zi| -wi) 

An + A2n,2n 



2(|2i|-u>i) (27) 



is an upper bound for the GM concurrence of X. Since 
2(|zi| — Wi) is also a lower bound for the concurrence, it 
is the exact value of the GM concurrence. 



Note that for the above decomposition to work we had 
to assume that |z;i| > w±. We now turn to the case 
\zi | < wi . We seek to show that all such density matrices 
are biseparable. We consider two different scenarios. 

(b) V ai&i > Wi- In this case the matrix X can be 
decomposed to matrices similar to the previous case. 



X 



i>l 



(28) 



where S[ is an Xu matrix whose two-qubit counterpart 
reads 



R' 



aiTi ziTi 
ai Zi 



z{Ti 



, Ti 



biT t 



Wl 



(29) 



Since |zi|T; < \fa^>i and \zi\ < \/aibiTi, we can invoke 
Lemma 1 to confirm R[ is biseparable for all i. The 
fact that is not normalized does not interfere with 
the proof of biseparability as one can always factor out 
Tr[fi-]. Now we focus on the last case. 

(c) V ai&i < wi. In this case we divide our matrix into 
two positive semi-definite matrices X = K\ + L\. 



I ait 



Ki = 



z x t\ 



a 2 r 



z 2 r 



z*r b n r 



bor 



\ztt 
( air 



(30) 



a 2 t 



z 2 t 



bit/ 
zir \ 



where 



\zir 



t = 



z* 2 t 



Wl 



wi + y/aibi 



CL n t Z n t 

z*t b n t 



b 2 t 

b x rj 

, and r = 1 — t. 



(31) 



Note that since wi < 3^/aibi, then | > t > r. One 



can show that 



t^Jaibi 



rwi, 



(32) 



which guarantees that Ki falls in the category of case 
(a). Since 



t\zi\ < rwi, 



r\zj | < (t — r) V aibi + rw 3 , 



(33) 



MM M 



K\ is biseparablc. Regarding matrix L%, since 
r|zi| < tWi, 

n 

< r\/ aibi + * £ y/^A, (34) 

it does not fall in the category of case (a) and thus belongs 
to either case (b) or case (c). If it falls in the category 
of case (b) then we can conclude that it is biseparablc. 
If not, we divide L\ into two matrices L± = K 2 + £21 
as before. Each time we divide a matrix in this way the 
trace of the remaining part is strictly smaller than the 
trace of the step before: Tr(L^) < 0.75*. Thus, we can 
write the matrix X as a convex sum of biseparable states 
and a remaining part that can be made arbitrarily close 
to zero. Therefore matrix X is a biseparablc matrix. 
This completes the proof for all X-matrices. Therefore, 
we have proved that the GM concurrence of a N-qubit 
X-matrix is 

C GM (X) = 2max{0,\z l \~w l }, l<i<n (35) 

IV. TV QUBITS IN TV CAVITIES 

The simple algebraic formula that was derived in the 
previous section opens up the possibility of analytically 
studying the entanglement dynamics of TV qubits in dif- 
ferent scenarios. The only restriction is that the density 
matrix of the system of interest remains an X-matrix. In 
this section we study the dynamics of the TV-partite en- 
tanglement between TV two-level atoms, each residing in 
its own single mode cavity, see Fig. [TJ The interaction 
between each atom and its partner field is described by 
a Jaynes-Cummings Hamiltonian. The Hamiltonian for 
the total system reads 

H <ot = £ y ^ + g(4a® + ai af) + ua\ ai , (36) 

i 

where loq , tu are the frequencies of the atoms and cavi- 
ties, respectively and g is the coupling constant between 
each atom and its partner cavity. For simplicity we have 
assumed that all the atoms and all of the cavities have 
the same resonance frequencies, ujq and w, respectively. 
We assume that all of the atoms are initially in an N- 
qubit GHZ state, sharing TV-partite entanglement, and 
the cavities are all in the vacuum state. 

We present a detailed analysis only for the three qubit 
case but the generalization to TV qubits is straightfor- 
ward and we present the result for TV qubits at the end. 
Our system is a generalization of the scenario previously 
studied by Yonag et al. [22|, so we follow the notation 
developed there. In fact, this generalization has already 
been employed in the study of multipartite entanglement 
in by Man et al. [25j , where the entanglement of differ- 
ent biseparations of the system were studied. Although 
separability of a particular bipartition of the system is a 



FIG. 1. A schematic representation of the setup where the 
remote atoms are interacting with their own cavities. There 
is no interaction between remote parties. The qubits are as- 
sumed to be initially maximally entangled and the fields are 
in individual separable vacuum states. 



sufficient condition for biseparability it is not a necessary 
condition. Thus a quantitative examination of TV-partite 
entanglement was missing; this is the analysis that we 
seek to provide using our formula. For simplicity from 
now on we refer to GM concurrence simply as concur- 
rence. 

We first consider an initial state which is an analog 
of the so-called |$,a) states for three atoms [22j . The 
initial state of the system reads 

= (cosa|e,e,e) + sin a\g, g,g}) <8> |0,0,0). (37) 

where |e) denotes the excited state of a two level atom 
cr z \e) = |e) and \g) denotes the ground state <J z \g) = 
— \g) . All three cavities are assumed to be in their ground 
states initially. Here we use a slightly different notation 
than we used in the previous section, where and 1 were 
used for bits of a qubit. From now on, and I will be 
used for the vacuum state and a singly excited state of 
the cavity. For each atom and cavity the two states |e, 0) 
and I g,0) will evolve according to 

U(t)\e,0)=A t \e,0)+B t \g,l), 

U(t)\g,0) = \g,0), (38) 

where U (t) is the Jaynes-Cummings propagator. If the 
atoms are in resonance with the cavities (ljq = w), then 
one can show that 

A t = e lu,t cos gt, B t = -ie~ wt sin gt. (39) 

Thus our initial state evolves to 

|#(t),a) =cosa {A t \e,0) + B t \g,l)) 
®(.A t |e,0>+B t |ff,l» 
8>(A t |e,0>+ B t |ff, 1» 
+ sina Iff, 0)® |ff, 0)8) 0). (40) 

By tracing out the cavities we find the density matrix of 



FIG. 2. Qjv(t) for different number of qubits. N — 2 (dashed 
line), N = 10 (dottedd line), N = 20 (dashed-dotted line), 
and N = 100 (solid line). The initial states are |$jv, j)- 



FIG. 3. i?jv(t) and T N (t) for different number of qubits. N = 
2 (dashed line), N = 10 (dottedd line), iV = 20 (dashed- 
dotted line), and TV = 100 (solid line). The initial states are 

l^,|>. 



the three atoms, 

/ ai 



/5<J> 3 (*) 



ft 2 



"2 



ft 2 



(41) 



6l/ 



where 



ai = cos 2 a|A t | 6 , 
b\= sin 2 a + cos 2 a|i?t| 6 , 
z\ = sin a cos aA 3 , 
ft2 = cos 2 a\A 2 B t \ 2 , 
6 2 = cos 2 a|^ tJ B 2 | 2 . 

The concurrence of this density matrix is 



(42) 



C 3 (t) = max{0,Q 3 (t)}, 

Q 3 {t) = 2 cos 2 a I cos 3 . gt I (| tana| - 3| sin 3 gi|) . (43) 

If, instead of the three qubit-cavity systems, we had N 
qubit-cavities then one can show that the concurrence of 
the N qubits would be 

C N (t) = max{0,Q w (i)}, 
Q N {t) = R N (t) T N (t), 
R N (t) = 2 cos 2 a\ cos N gt\, 

T N (t) = I tana| - - 1)| sin^ gt\. (44) 

For N = 2, we recover the formula given by Yonag et al. 
. We plot Q N , T N and R N in Fig. [2] and Fig. [3] for 
different numbers of qubits ranging from 2 to 100. As 
N increases, the positive part of Tn approaches a step 
function (when T/v < 0, the concurence is zero), whereas 
Rn remains non- negative and only vanishes momentar- 
ily. Thus TV determines when entanglement vanishes 



completely and Rn determines the magnitude of concur- 
rence. Notably = Rn-i\ cos gt\, which guarantees 
that a as N increases the bulk of the initial concurrence 
decays faster. The time it takes for the concurrence to 
decay by a factor of two can be approximated by 



9Ti 



arccos ■ 



V2' 



(45) 



The other time of interest is when the concurrence van- 
ishes completely. For N > 3, we can directly confirm the 
phenomenon of multipartite entanglement sudden death 
(ESD). To find the time of ESD, T ESD , we solve the 
equation Qjv = which gives 



Tj 



ESD 



1 

— arcsm 
9 



tana 



2 N-i _ 1 



lim Tesd = t^- 
N->oo bg 



(46) 



In the limit of large N, the time of ESD is independent 
of the initial concurrence; it depends only on the com- 
mon coupling strength g. In Fig. [4] we have plotted the 
ESD time against the number of qubits. When the initial 
amount of concurrence is considerable the ESD time be- 
comes smaller as N increases. However, when the initial 
concurrence is small the Tesd increases and converges to 
the limit Tesd = 7r/6g as TV tends to infinity. It is no- 
table that even for a very large N, although the bulk on 
initial entanglement decays very rapidly, the state of the 
two-level atoms remains entangled and becomes bisepa- 
rable only after Tesd- 

We also consider an initial state which is the three- 
qubit analog of the so-called |\&) state. 

I*, a) = (cosa|e,e,5) + sin a\g , g , e)) ® |0,0,0) (47) 



This state evolves according to 



(48) 



|*(t),a) =cosa (A t \e,0) + B t \g,l)) 

®(A t \e,0) + Bt\g,l))®\g,0) 

+ sma 0)8) |s, 0) <8> (A t \e, 0) + B t \g, 1)). 
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FIG. 4. ESD time versus number of qubits for different initial 
states |$jv, a). From the bottom up tana = 0.01, 0.1, 0.2, 0.5, 
and 1 respectively. 



FIG. 5. Qjv(t) versus time for different numbers of qubits. 
N = 2 (dashed line), N = 10 (dottedd line), iV = 20 (dashed- 
dotted line), and N = 100 (solid line). The initial states are 



The reduced density matrix of the three atoms in the 
basis {|e,e,e), \e,e,g), \e,g,e), . -.,\g,g,g)} is 



(*) = 



\ 



"2 



?2 



«4 



-2 



(49) 



where 



a-2 
z-i 
a 4 

l>2 

bi 



cos 2 a\A t \ 4 

sin a cos <xA\A* t 

c(x?a\A t B t \ 2 

sin 2 a\A t \ 2 

cos 2 a\A t B t \ 2 

sin 2 a|B t | 2 +cos 2 a|B t 



The concurrence is given by 

Cz(t) = 2| sin a cos a cos 3 | . 



(50) 



(51) 



This state does not experience ESD due to the fact that 
(e,e,e\py\e,e,e), (e,g,e\py\e,g,e), and (g, e, e\p^\g, e, e) 
remain zero for all times. To show this, we note that our 
initial state is a superposition of two possibilities. Either 
two of the atoms, {1,2}, are excited and the other atom, 
{3}, is in vacuum state, or that single atom is excited 
and the atoms {1,2} are in their ground states. Since all 
cavities are initially in their ground states, if an atom is 
in its ground state initially it will always remain there. 
But the three diagonal terms that we referred to require 
the atom {3} and at least one of the other two atoms 
be excited simultaneously. Since this is forbidden, all 
of these matrix elements remain zero. Thus V a i^i = 
y/a^b 3 = Vfl4^4 = 0, and C3 = 2|z 2 |. This argument can 
be generalized for TV > 4 to all of the GHZ states except 
for |$jv,a), because in the |4>jv,a) state all of the atoms 



can be initially excited. Therefore none of the GHZ states 
except \<f>N,a) experience ESD. The concurrence of all 
such GHZ states is given by 



Cjy(i) = 2| sin a cos a cos gt\ 



(52) 



In Fig. [5] we plot Qn for different numbers of qubits 
ranging from 2 to 100, The initial state is I^jv, f )• As 
expected from the concurrence formula, ESD does not 
happen for any N and only vanishes momentarily. Still, 
as N increases the width of the periods of considerable 
concurrence decreases. 



V. CONCLUSION 

We have found an algebraic formula for the genuinely 
multipartite (GM) concurrence of N-qubit density ma- 
trices that can be written as X-matrices in an orthonor- 
mal product basis. This development allows iV-partite 
entanglement to be quantified for such states. The for- 
mula opens up the possibility of studying entanglement 
dynamics of ./V-qubit states in different scenarios, as long 
as the X-form of the density matrix is preserved. This in- 
cludes scenarios such as local decoherence channels that 
will be studied elsewhere. 

Using the concurrence formula, we have studied the dy- 
namics of iV-partite entanglement of N two-level atoms 
interacting with partner single-mode cavities. We showed 
that the so-called |$j\r, a) states experience ESD whereas 
the I 'J at, a) states do not experience ESD. Algebraic for- 
mulas for the concurrence of the two different case were 
presented. Since the entanglement dynamics is periodic 
we focused on how fast the initial entanglement vanishes 
in the initial cycle before reviving to its initial value at 
gt = 2ir. It is observed that although the bulk of initial 
concurrence decays faster as the N increases, the same 
conclusion can not be drawn about the exact moment at 
which the entanglement disappears, Tesd- In fact for 
the states that have small initial entanglement, Tesd is 
prolonged as the number of qubits increases. 
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